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ABSTRACT 

This  research  indicates  that  a  knowledge  of  angle  of  attack  (a) 
and/or  gust  perturbation  angle  of  attack  (a  )  is  of  prime  importance 
for  a  longitudinal  autopilot  designed  to  alleviate  the  effect  of 
vertical  gusts.  However,  accurate  measurement  of  a  and  a  is  extremely 
difficult.  This  research  indicates  that  the  readily  measurable  normal 
acceleration  factor  (n)  and  pitch  rate  (q)  can  be  used  to  provide  a 
and  a  information  indirectly. 
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SYMBOLS 


a         =  Linear  acceleration  of  aircraft  center  of  mass  in 
cz         z  direction,  ft/sec. 

c         =  Length  of  mean  aerodynamic  chord,  ft. 


c 


8c 
m 


m  dr\ 

n  o 


CL 


=  3CL 


a  Sa  Q 

g  =  Local  acceleration  of  gravity,  ft/sec. 

h  =  Altitude,  (ft). 

e2  =  Mean-square  value  of  output  of  system. 

K1->K4  =  Control  system  gains. 

L  =  Gust  characteristic  length,  ft. 

m  =  Aircraft  mass,  Slugs. 

— ^  =  a  /g,   normal  acceleration  factor,  ft/sec 2/g. 
z 

&  m  ^unallev.  "  &   '  ("/secVg)2. 

unallev.  =  Value  of  n2  with  no  autopilot  utilized,  ( ft/sec 2/g)  2. 

q  =  Pitch  rate,  rad/airsec. 

q_  =  Pitch  rate  due  to  vertical  gust,  rad/airsec. 

U  =  Aircraft  forward  velocity,  ft/sec. 

u  =  Aircraft  forward  velocity  perturbation,  ft/sec. 

xyz  =  Stability  axis  system. 

w  =  Velocity  perturbation  along  z  axis,  ft/sec. 

w  =  Gust  vertical  velocity,  ft/sec. 

a  =  w _/U  ,  angle  of  attack  perturbation  due  to  gust, 
measured  at  aircraft  center  of  mass,  rad. 


n  =  Elevator  deflection,  measured  from  trim,  rad. 

n2  =  Mean-square  value  of  elevator  deflection  (Rad)  2 

6f  =  DLC  (direct  lift  control)  deflection. 

6.  =  Throttle  deflection. 

Q  =  Spatial  frequency,  rad/ft. 


w 


ftU  t*  =  ftU  5ff  i  Frequency,  Rad/AirSec. 
o 


o  =  Root  mean-square  value  of  vertical  gust  velocity,  ft/sec, 

6         =  Pitch  angle  perturbation,  Rad. 
p         =  Atmospheric  density,  Slugs/Ft3. 
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I.   INTRODUCTION 

A.   LONGITUDINAL  AUTOPILOTS 

1.  Basic  Ideas  in  Longitudinal  Control 

Since  the  beginning  of  powered  flight,  methods  have  been  sought 
to  reduce  pilot  effort  in  controlling  aircraft  and  to  automate  various 
aspects  of  aircraft  control.  The  autopilots  and  stability  augmentation 
systems  of  today  are  realizations  of  such  research. 

Autopilots  and  stability  augmentation  systems  are  usually 
classified  as  longitudinal  or  lateral.  As  the  names  imply,  longitudinal 
(lateral)  autopilots  are  concerned  with  the  longitudinal  (lateral)  flight 
regulation.  This  dual  classification  is  justified  in  flight  conditions 
where  the  existence  of  pure  longitudinal/lateral  motion  is  ensured. 

Blakelock  [Ref .  1]  offers  a  very  readable  introduction  to  auto- 
pilot design.  He  discusses  three  basic  types  of  longitudinal  controllers, 
each  designed  for  different  longitudinal  operations:  a  displacement 
autopilot  where  aircraft  pitch  angle  is  fed  back  and  compared  with  a 
reference  pitch  angle;  a  pitch  orientational  autopilot  where  aircraft 
pitch  rate  is  fed  back  and  compared  with  some  desired  pitch  rate;  and  an 
acceleration  control  system  where  the  acceleration  of  the  center  of 
gravity  is  fed  back  and  compared  with  a  reference  acceleration.  A  good 
treatment  of  the  feedbacks  that  may  be  used  in  longitudinal  control 
systems  can  be  found  in  Ref.  2,  Sec.  II.  Table  I  (obtained  from  Ref.  2) 
indicates  the  different  feedbacks  which  can  be  utilized  to  control 
longitudinal  motion. 


The  research  summarized  here  is  concerned  with  autopilots  designed 
to  alleviate  the  effect  of  vertical  gusts  on  a  large  jet  transport  in 
steady,  wings-level  flight.  Two  flight  conditions  are  considered,  high 
altitude  cruise  and  landing  approach.  By  reducing  the  aircraft's  gust 
response  it  is  possible  to  provide  for  a  smoother  ride,  aircraft 
structural  fatigue  problems  can  be  alleviated,  and  the  ability  of  the 
aircraft  to  function  as  a  gun  platform  can  be  imporved. 
2.  Longitudinal  Control  to  Alleviate  Gusts 

Appendix  B  offers  an  introduction  to  the  statistical  design 
concepts  utilized  in  this  study.  As  Table  I  indicates  many  feedback 
parameters  such  as  pitch,  pitch  rate,  center  of  gravity  acceleration, 
etc.,  can  be  utilized  as  feedbacks  in  a  controller  designed  to  reduce 
gust  response.  Hess  [Ref .  3]  discussed  a  two-parameter  gust  alleviation 
system  shown  in  Fig.  1.  He  analyzed  this  system  using  a  performance 
index  given  by: 


Zti  *-*      ag   'a    '.~    ^  ag   'a   _  •* 

oo      =       g   S=JCl)      oo     =       g   s=jw 


where  w  is  the  non-dimensionalized  frequency  and  the  transfer  functions 

—  (s)  and  —  (s)  are  obtained  from  the  aircraft's  equation  of  motion  as 
«g       ag 

discussed  in  Appendix  C.  $   (w)  is  the  non-dimensionalized  power  spectral 

ag 

density  of  the  gust  field  and  was  obtained  from  Equation  1  found  in 
Appendix  B.  Hess's  study  utilized  a  large,  rigid  jet  transport.  The 
characteristics  for  the  two  flight  conditions  considered  may  be  found  in 
Appendix  A.  His  results  are  shown  in  Table  II. 
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Hess  also  proposed  an  optimal  control  system,  Fig.  2,  which 
showed  considerable  improvement  over  the  two-parameter  feedback  system 
above.  For  analysis  of  this  system  Hess  used  a  performance  index  given 
by: 

I  =  £^  S  [c^2(t)  +^2(t)]dt 

o 

The  aircraft  was  the  same  as  used  in  analysis  of  the  two-parameter  system. 
Table  III  shows  results  of  his  study  for  a  q/r  equal  to  10,  a  value  was 
selected  by  a  root-square  locus  technique.  Note  the  improvement  in 


An2  which  Hess '  s  optical  system  maintains  over  the  two- 

n  unalleviated 

parameter  system  of  Table  II. 


B.  BASIS  OF  THIS  RESEARCH 

A  major  obstacle  to  the  implementation  of  Hess's  optimal  gust 
alleviation  system  is  the  difficulty  in  measuring  a  .  Attempts  to  design 
instruments  that  will  give  satisfactory  on-line  measurements  of  a  have 
been  unsuccessful.  This  report  investigated  the  effect  of  eliminating 
a     information  on  the  performance  of  the  alleviation  system.  Figure  3 
shows  Hess '  s  optimal  gust  alleviation  system  with  gust  information 
removed.  This  is  the  system  utilized  in  this  study. 
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II.      SOLUTION  PROCEDURES 

A.  METHODS  OF  SOLUTION 

Optimal  feedback  design  can  be  classified  as  free  configuration  or 
fixed  configuration.  The  latter  is  often  referred  to  as  parameter  opti- 
mization. In  the  free  configuration  procedure  the  designer  is  not 
restricted  in  any  way.  The  final  design  depends  only  upon  the  input  and 
the  form  of  the  performance  criteria.  With  the  parameter  optimization 
techniques  the  configuration  is  chosen  using  feedbacks  which  are  readily 
obtainable  and  then  the  feedback  parameters  are  optimized  to  rrdnimize 
the  performance  criteria. 

As  noted  in  Ref .  2,  for  autopilot  design  the  parameter  optimization 
technique  is  probably  the  superior'  method  considering  current  state  of 
the  art.  This  technique  yields  an  autopilot  which  fits  the  "real  world" 
situation.  This  research  utilized  the  parameter  optimization  procedure 
since  the  form  of  the  longitudinal  controller  was  prespecified  (Fig.  3) . 

B.  PERFORMANCE  INDEX 

The  performance  index  chosen  for  the  optimization  procedure  was 
given  by: 


2tr 


_c°  °°  ri 

£$n(w)du)  +  £  $  (co)dwl  =  n2+ n2  (1) 


where  u  denotes  the  non-dimensionalized  frequency  of  the  random  gusts. 
This  performance  index  is  identical  to  that  used  by  Hess  in  his  evaluation 
of  the  optimal  controller  with  q/r  =  10.  It  is  a  sum  of  the  mean  square 
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normal  acceleration  factor  and  the  mean  square  elevator  deflection.  The 
latter  was  included  to  prevent  large  elevator  excursions.  Equation  (1) 
can  be  manipulated  as  indicated  in  Appendix  B  to  a  more  usable  form: 


1  -  h  S  V;>  l=-0->  |2d"  +  ST  S*^M  I2-(W  \2^         (2) 

°°       3     g  °°  3     g 


$   is  obtained  from  Equation  (1)  in  Appendix  B  by  substituting  to  =  ftu  t*, 
This  non-dimensionalized  form  is: 


$  C)   =   JL£0   , _i (3) 

agV  '    LwUq2    a2   +    (C-2/4L2) 

The  transfer  functions  —  (jw)  and  —  (jw)  are  obtained  from  the  air- 

g      g 

craft's  longitudinal  equations  of  motion,  the  equation  for  normal 
acceleration  and  the  mathematical  equation  of  the  autopilot.  Appendix  C 
shows  how  this  was  done. 

For  the  purposes  of  determining  the  aerodynamic  gust  effects,  the 
aircraft  was  idealized  as  a  segment  of  the  x-axis.  This  idealization 
has  been  found  to  give  good  results  as  long  as  the  shortest  vertical- 
gust  wave  length  is  greater  or  equal  to  eight  times  the  aircraft  length 
[Ref .  4,  pg.  323] . 

C.   INTEGRATION  OF  THE  PERFORMANCE  INDEX 

Substitution  of  Equation  (3)  and  the  transfer  functions  into 
Equation  (2)  produced  an  equation  whose  form  is  frequently  found  in  system 
analysis.  Tables  for  solution  of  these  equations  are  readily  available, 
e.g.  Ref.  5. 
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Appendix  D  shows  how  a  computer  program  was  written  to  solve  for 
the  performance  index,  which  in  this  case  took  on  the  functional  form 
of  I/(KlfK2,K3,L  )  . 

D.   OPTIMIZATION 

Integration  of  the  performance  index  for  a  given  characteristic  gust 
length  (L  )  produced  a  function  which  contained  the  variables  Kl,  K2, 
and  K3,  i.e.  the  feedback  gains.  The  function  obtained  could  be  thought 
of  as  a  performance  index  surface  with  dimensions  Kl,  K2  and  K3.  The 
optimum  solution  would  be  found  where  the  surface  produced  a  minimum 
value.  A  computer  program  was  written  to  solve  for  this  minimum  by  a 
method  of  steepest  descent.  Appendix  D  indicates  how  this  was  accomplished, 
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III.   RESULTS  AND  DISCUSSION 

A.   PERFORMANCE  COMPARISONS 

The  results  of  the  computer  computation  are  presented  in  Table  IV. 
The  performance  surface  at  most  points  was  very  shallow  making  it  diffi- 
cult to  rapidly  find  a  minimum  value  on  the  surface.  It  was  found  that 
the  performance  index  was  a  weak  function  of  K2,  but  a  strong  function 
of  Kl,  K3  and  L  .  Kl,  K2  and  K3  were  found  to  be  very  dependent  on  the 
flight  condition,  i.e.  cruise  or  landing  approach. 

Comparing  Table  IV  with  Tables  II  and  III  shows  that  the  three- 
parameter  system  has  the  largest  performance  index.  Figures  5  and  6 
indicate  that  the  An^/n2  performance  initially  compares 

favorably  with  the  optimal  and  the  two-parameter  system,  but  the  per- 
formance rapidly  deteriorates  as  L  increases.  This  indicates  that  a 
knowledge  of  a  becomes  increasingly  important  as  L  gets  larger. 

Note  that  in  the  two-parameter  system  of  Fig.  1  the  feedbacks  n  and 
q  give  a  knowledge  of  a  and  a  .  This  can  be  seen  from  the  equation  for 
the  normal  acceleration  factor  which  can  be  written  as: 


D(a  +  a  )  =  q  -  |p>  n 

3         o 


B.   TRANSIENT  PERFORMANCE 

Figures  6  and  7  present  the  locus  of  roots  of  the  characteristic 
equation  obtained  when  the  feedback  gains  Kl  and  K3  are  set  at  values 
obtained  for  a  gust  characteristic  length  of  1000  feet  while  K2  was 
allowed  to  vary.  This  could  be  done  because  the  performance  index  was 


15 


a  weak  function  of  K2.  A  value  of  the  variable  K2  could  be  selected 
which  would  give  the  three-parameter  system  better  transient  performance 
than  the  two-parameter  system,  but  this  transient  performance  was  still 
inferior  to  the  optimal  control  system.  A  new  root  locus  diagram  would 
have  to  be  drawn  for  each  gust  length  and  flight  condition.  In  an  effort 
to  circumvent  the  drawing  of  multiple  root  locus  diagrams,  the  feedback 

gains  Kl  and  K3  were  set  at  values  obtained  for  L  equal  to  3000  feet. 

*  w 

It  is  found,  however,  that  by  doing  this  performance  was  unacceptably 
degraded,  particularly  in  the  lower  ranges  of  L  . 

C.  STABILITY 

It  was  found  that  the  three-parameter  system  was  stable  throughout 
the  range  of  feedback  gains  considered  here. 

D.  POSSIBLE  ALTERNATE  SOLUTION 

Comparing  the  three-parameter  longitudinal  controller  of  this 
research  with  the  optimal  controller  (Fig.  2)  and  the  two-parameter 
system  (Fig.  1)  it  appears  that  a  is  an  important  parameter  for  a  gust 
alleviation  system.  Since,  using  present  technology,  it  is  not  possible 
to  measure  a  aboard  the  aircraft  some  alternative  method  of  obtaining 
this  information  must  be  employed. 

It  was  noted  previously  (Sec.  Ill,  Part  A)  that  a  knowledge  of  n  and 
q  would  provide  information  on  a  and  a  .  This  would  explain  why  the  RMS 
performance  of  the  two  parameter  feed  back  system  tends  to  approach  the 
RMS  performance  of  the  optimal  controller  as  the  characteristic  gust 
length  increases. 
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It  is  felt  that  by  cxmbining  the  optimal  controller  and  the  two- 
parameter  system  a  gust  alleviation  system  could  be  built  that  would  have 
performance  comparable  to  the  optimal  controller.  Indeed,  Hess  [Ref .  6] 
has  shown  this  to  be  the  case.  The  modified  gust  alleviation  system 
appears  in  Fig.  4.  This  system  would  have  the  advantage  of  utilizing 
feedback  information  that  is  readily  obtainable  on-board  the  aircraft. 
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IV.  CONCLUSIONS 

1.  The  longitudinal  gust  alleviation  system  evaluated  by  this  report  is 
a  stable,  usable  system  and  could  be  designed  to  provide  better  transient 
performance  than  the  two-parameter  system  of  Fig.  1.  However,  the  two- 
parameter  system  and  the  optimal  controller  of  Fig.  2  both  exhibit  better 
RMS  performance. 

2.  a  provides  important  information  for  a  gust  alleviation  system. 
Knowledge  of  a  becomes  increasingly  important  as  the  gust  characteristic 
length  increases,  i.e.  as  the  random  gust  field  departs  from  white  noise. 

3.  lb  provide  for  a  gust  alleviation  system  which  uses  readily  obtainable 
information,  and  still  produce  satisfactory  performance,  the  system 
presented  in  Fig.  4  is  proposed. 
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Table  II 


TWO-PARAMt'l 

'LR  LONGITUDINAL  AUTOPILOT 

Cruising 
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An2 

Kl 

Pt/Airs2/g 

K2 
Rad/Rad/Airs 

L 
w 

;   Ft 

(Ft/Sec2/g)2 

unallev , 
(Rad)2 

0 

w 
(Ft/Sec) 

.36 
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500 

.0466 

.00282   .27 

10 

.36 

475 

1000 
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.00153   .37 

10 

.36 

400 
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.0151 

.000818   .41 

10 

.36 

375 

3000 

.0139 

.000501  .24 

10 

.36 

350 

4000 
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.000435   .44 

10 

.36 

350 

5000 

.00637 

.000348   .44 

10 

.36 

350 

6000 
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.000290   .44 

10 

Landing  Approach 


Kl 
.74 
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.74 
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.74 
.74 
.74 
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K2 

L 
w 

n2 

n2     unallev. 

o 
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500 

.0384 

.0138 

.21 

10 
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1000 

.0209 

.00690 
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10 

175 
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.0109 
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10 

175 

3000 
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10 
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4000 
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10 
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Table  IV 


^ 

Cruising  Flight 

LJ   \_-Wl  U^  \is-/-t_U-lLJ 

Kl 

K2 

K3 

"n^ 

n2 

to2" 

Ft. 

Rad/Rad 

Rad/Rad 

Rad/Rad 

(Ft/! 

(Rad)2 

n2 

Airs 

UJlcL 

500 

1.59 

688.0 

-2.57 

.0461 

.0001 

.276 

1000 

1.60 

688.0 

-2.57 

.0324 

.0001 

.259 

2000 

0.795 

688.0 

-3.08 

.0207 

* 

.201 

3000 

0.568 

688.0 

-3.07 

.0150 

* 

.176 

4000 

0.452 

688.0 

-2.99 

.0120 

* 

.143 

5000 

0.393 

688.0 

-2.96 

.0100 

* 

.123 

6000 

0.447 

688.0 

-2.99 

.0085 

* 

.115 

a   =  10  Ft/Sec 

*  -  Value  less  than  .0001 


Landing  Approach 


h, 

Kl 

K2 

K3 

n2 

n2 

An5 

Ft. 

Rad/Rad 

Rad/Rad 
Airs 

Rad/Rad 

(Ft/Sec2/g) 2 

(Rad)2 

n2 
un« 

500 

0.651 

400.0 

-1.00 

.0356 

.0011 

.262 

1000 

0.785 

400.0 

-1.13 

.0219 

.0008 

.270 

2000 

0.255 

400.0 

-0.996 

.0124 

.0004 

.250 

3000 

0.165 

400.0 

-0.939 

.0086 

.0003 

.265 

4000 

0.165 

400.0 

-0.939 

.0066 

.0003 

.236 

5000 

0.098 

400.0 

-0.974 

.0054 

.0002 

.250 

6000 

0.085 

400.0 

-0.953 

.0045 

.0002 

.250 

o   =  10  Ft/Sec 
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Two  Parameter  Gust  Alleviation  System 
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FIGURE  2 
Optimum  Gust  Alleviation  System 
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FIGURE  3 
Three  Parameter  Gust  Alleviation  System 


FIGURE  4 
Proposed  Gust  Alleviation  System 
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FIGURE  7 
Three  Parameter  Gust  Alleviation  System 
ROOT  LOCUS 
CRUISING  FLIGHT 
Kl  =  .785,  K3  =  -1.13 
1^  =  1000  ft. 
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FIGURE  8 
Three  Parameter  Gust  Alleviation  System 
FOOT  LOCUS 
LANDING  APPROACH 
Kl  =  .785,  K3  =  -1.13 
L^  =  1000  ft. 
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APPENDIX  A 
Aircraft  Characteristics 

The  aircraft  utilized  in  this  report  was  the  large,  rigid,  jet 
transport  used  by  Etkin  to  Illustrate  application  of  theory  he  presented 
in  Chapter  6  of  Ref .  4.  The  cruising  flight  data  were  obtained  from 
Ref .  4  and  the  landing  approach  data  were  obtained  frcm  Ref.  7. 


Parameter 


Cruising  Flight 


Landing  Approach 


W 

S 

A 

Alt. 

U 
o 

P 

c/2 

y  =  m/psc/2 

t*  =  c/2U 
'  o 

ig  =  B/psc/2 


°L 
c 


m 


m- 


1000,000  lbs. 

1667  sq.  ft. 

7 

30,000  ft. 

500  mph.  or  733  fps. 

0.000889  slugs/ft.3 

7.7  ft. 

272 

0.0105  sed. 

1900 

0.25 

0.0188 

-4.9 

0 

-0.488 

-4.2 


1000,000  lbs. 

1667  sq.  ft. 

7 

0  ft. 

200  mph.  or  294  fps, 

0.002377  slugs/ft.3 

7.7  ft. 

102 

0.0262  sec. 

710 

0.59 

0.0295 

-4.8 

0 

-0.478 

-4.2 
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Aircraft  Characteristics  (Continued) 
Parameter Cruising  Flight landing  Approach 

C  -22.9  -22.9 

m 

q 
c  oo 

z 

q 

C  -0.24  -0.24 

z 

n 

C  -0.72  -0.72 

m 


n 

c  . 

m 
n 


o  o 
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APPENDIX  B 
Randan  Gust  Fields  and  Their  Use  in  Autopilot  Design 

Aircraft  controllers  designed  to  behave  well  when  disturbed  by 
sinusoidal,  step  or  impulsive  gust  disturbances  often  perform  poorly 
when  subjected  to  actual  atmospheric  conditions.  This  has  led  to 
seeking  better  methods  of  defining  the  turbulence  or  gust  fields.  It 
has  been  discovered  that  if  the  atmospheric  disturbance  was  treated  as 
a  random  process  and  the  aircraft  controller  designed  using  statistical 
concepts  many  of  the  above  difficulties  could  be  eliminated.  Controllers 
could  then  be  analytically  designed,  and  built,  which  performed  well 
when  subjected  to  actual  atmospheric  conditions. 

Statistical  design  techniques  brought  with  them  many  new  terms  and 
concepts  not  used  in  the  era  of  the  so-called  classical  autopilot  design. 
Terms  such  as  stationary,  homogeneous,  isotropic,  frozen  turbulence, 
mean-square  error,  correlation  functions  and  power  spectral  densities 
are  now  included  in  the  analytical  designer's  vocabulary.  This  writer 
will  not  attempt  to  provide  any  mathematical  derivations  of  these 
concepts  and  will  not  even  discuss  them  at  length.  The  interested 
reader  is  referred  to  the  bibliography  found  in  Ref.  4,  pp.  339-340. 
A  very  readable  review  of  this  subject  is  found  in  Ref.  1,  Chapter  9. 

In  order  to  utilize  statistical  techniques  in  gust  alleviation 
design  an  expression  for  the  power  spectral  density  of  the  atmospheric 
turbulence  is  needed.  Such  expressions  have  been  developed  which  show 
good  agreement  with  spectra  measured  in  cumulus  clouds,  thunderstorms 
and  clear  air.  The  spectra  used  for  the  random  vertical  gust  field  of 
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this  report  is  found  in  Ref .  3,  given  by: 


VQ)  =  VLw  r+Srop  (1) 


where  a  2  is  the  mean  square  value  of  the  vertical  gust  velocity.  L  is 
the  characteristic  gust  length,  and  ft  is  the  spatial  frequency. 

Equation  (1)  can  be  written  in  non-dimensional  form  as: 

2 

wc         i 

*wg(w)  =  IT"  (52  +  (cV4L^)  (2) 


Although  this  expression  is  not  as  comprehensive  as  many  expressions 

found  for  <J>Wcf,  it  is  used  because  it  is  easy  to  manipulate  and  it  agrees 

fairly  well  with  measured  atmospheric  power  spectra.  A  form  of 

Equation  (1)  was  originally  suggested  in  Ref.  8.  a     has  been  found  to 

w 

vary  from  0  to  40  feet  per  second  and  L  from  500  to  6000  feet. 

w 

The  power  spectral  density  of  any  of  the  aircraft  longitudinal 
motion  variables  such  as  angle  of  attack  (a)  can  be  obtained: 


*„,>>  =*,»>>  I  £-Ml2  (3) 

act     wq   '  w    ' 

y    g 

Here  —  (jw)  represents  an  aerodynamic  transfer  function,  with  s  =  jw, 

W 

g 

and  relates  angle  of  attack  and  gust  vertical  velocity. 
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In  addition,  mean  square  values  of  longitudinal  output  variables 
can  be  obtained  as: 

CO 

cF=  ~-  C*   (a))du  (4) 


=  hS  w^i^i2*;         (5> 


g 


Equations  (3)  and  (5)  provide  the  necessary  tools  for  simple  analysis 
and  parameter  optimization  studies. 
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APPENDIX  C 
Determination  of  Transfer  Functions 


|  (»  and  £  (jw) 

g      g 


I.  longitudinal  Equations  of  Motion. 

The  aircraft  longitudinal  equations  below  were  obtained  from  Ref .  4, 
page  133: 

A. 

X  force:  (2yD  -  2CT  tane  -C  )u  -  C  a  +  C_  8  =0 
K      L     ox      x     L 


o         u 


Z  force:  (2CT   -  C  )u  +  (2uD  -  Cz  .D  -  Cz  )  a  -  [  (2y  +  C  )D 
o    Zu  a      a  z 

-  C^  taneQ]e  -  CZnn  =  0 


Moment: 


-q-  u  -  (C^.D  +  C^  )a  +  (i  D2  -  C^  D)9 

u       a      a       °       a 


q  -  D6  =  0        t*  =  l/u  I   =  c/2 

Substituting  a=a+a  /  q  =  q  +  q?  and  p  =  p  =  p  ; 

Neglecting  changes  in  forward  velocity  and  considering  the  X-force 
equation  to  be  identically  satisfied  yields  Equations  (1)  and  (2)  below. 
The  equations  are  shown  as  they  appear  after  the  Laplace  transform  has 
been  taken  with  zero  initial  condition  assumed: 


(2yS  -  C  s  -  C  )a  -  (2p  +  C  )q  -  C  H  =  [C  s  +  C   -  C  s]  a    (1) 
z.     z  zzn      z,     z     zq 

a      a  q  a      a     q 


tCm.s  +  Cm  )a  +   (iBs  "  Cm  )(5  -   [Cm-s  +  Cm  >»  (2) 

a  a  q  n  n 


[C     s  +  C        -  C     s]a 

m.  m  in         g 

a  a  q        3 
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note:   (a)  q  +  q  =  D6 

.   _  d(Mg)dt  _   uo 
%3   "  dt  dx    dt 


q  =  -Da 

(b)  The  a  terms  above  are  included  on  the  right  side  of  Equations 
(1)  and  (2)  to  indicate  their  use  as  forcing  functions. 

II.  Acceleration  Equation. 

The  equation  for  the  linear  acceleration  of  the  aircraft  center  of 

mass  in  the  z  direction  is  obtained  from: 

a   =  ft  +  pv  -  QU 
c 
z 

by  neglecting  velocities  in  the  y  direction,  assuming  small  disturbances, 

and  using  w  =  u  sina  Equation  (3)  is  obtained: 
o 

ac  =  J>  (Da  -  q)  (3) 


Now  one  can  define  the  normal  acceleration  factor: 

c    2u  *     A 

n(t)  =  -r  =  ST"  (q  -  Da) 


or  in  Laplace  form 


2U  2 
n(s)  =  ^-  I*   -  Sa>  (4) 
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III.  Mathematical  Model  of  Longitudinal  Control  System. 

The  equation  for  the  longitudinal  control  system  was  obtained  from 
the  block  diagram  of  the  system  and  is  given  by: 

Kla  +  K2q  +  (K3  -  TS  -  1) n  =  0  (5) 

Note  that  the  elevator-serro  combination  is  approximated  by  a  first 
order  system  with  a  time  constant  (T)  of  1/10  second.  T  in  the 
non-dimensionalized  form  is  given  by: 


2TO 
T*  =  T/t*  =  — • 


IV.  Simultaneous  Solution. 

From  Equations  (1) ,  (2) ,  (4) ,  and  (5)  a  system  of  equations  is 
obtained  of  the  form: 

Ana  +  A12q  +  A^n  +  A^n  =  b^ 

A21a  +  A22q  +  A^n  +  A^n  =  b^  (6) 

A^a  +  A32q  +  A^n  4-  A34n  =  b3ag 

A4la  +  A42^  +  A43n  +  A44n  =  b4ag 


where 


A   =  2yS  -  cz  s  -  Cz        A   -  Cfo.s  -  On 
J.J-  a      a         ZJ-     a      c 

A12  "  -2v  -  \  *22  "  V  "  \ 

A13  "  "  \ 

A14  =  ° 

b,=Cs  +  C   -  C  s 
1    z.     z     z 

a      a     q 


Ao^  ~  "On-   ~  On 

A24  =  ° 

bn  =  C     s  +  c 

2        m.           m 
a              a 

C     s 

m 

q 
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*3L 

=  Kl 

^32 

=  K2 

^3 

=  K3 

A34 

=  0 

b3  = 

0 

T*s 


A41  ="2U02s/gc 


A42 

=  2U 
o 

Vgc 

A43 

=  0 

A44 

=  -1 

b4  = 

'  0 

The  transfer  function  n(s)/a  (s)    is  given  by: 

a     '  '. 


n 


.(S) 


s=jw 


All 

^2 

A13 

bl 

A21 

h.2 

^3 

b2 

hi 

^2 

A33 

b3 

A41 

A42 

A43 

b4 

All 

A12 

A13 

A14 

^1 

A22 

A23 

A24 

A31 

A32 

A33 

^4 

A 
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M2 


A 


43 


A 


44 


(7) 


n(s)/a  (s) 
y 


is  obtained  in  a  similar  manner. 


S=JU 
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APPENDIX  D 
Computer  Program 

I.  The  computer  program  was  written  to  provide  a  solution  to  a  performance 
index  given  by: 


CO  CO 

n   2tt  sJ   ag'a     '      2tt  — >  a„  'a  J   ' 


(1) 


where 


I   (w) 
ag 


a2c 

ru2 

w  o 


1 

W2  +  C^ 


4L^ 
w 


(2) 


H  (Jw)  |2  =  xl(j(o)2  +  x2(jco)  +  x3 

x4(joj)3  +  x5(jw)2  +  x6(jw)  +  x7 


(3) 


n  ,.\ 
— (3«) 

g 


Zl(joj)3  +  X2(jgj)2  +  z3(joj)  +  z4 


x4  (ju) 3  +  x5  ( jw)  2  +  x6  (ju)  +  x7 


(4) 


In  Equations  (3)  and  (4) ,  xl  through  x7  and  zl  through  z4  contain 
the  variables  ¥1,  Y2,   and  K3  and  coefficients  of  these  variables.  The 
coefficients  are  functions  of  the  aircraft  stability  derivatives  and 
the  gust  characteristic  length. 

It  can  be  seen  that  the  performance  index  can  be  written  in  the  form: 


co  ^ 

i  BI  r  c(joj)c(-jo))  dw 

n   2tt  J  d(jw)d(-jw) 


(5) 
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This  integral  is  evaluated  for  n  =  1  through  n  =  10  in  Ref .  5.  For  the 
performance  index  of  this  problem  1=1.  where 

c|  (-dfe-Kipd^)  +  (cl-a^Cg)  d0dl(34+  (c|-2cQc2)  d0a3d4+cg  (-d1d|+d2d3d4 ) 
I  —  -     (6) 

2d0d4  t-d^l^d^d^) 

A.  *.  *»  A 

c(ju>)  +  c3(jw)3  +  c2(joj)2  +  c^jw)  +  Cq 

d(jw)  =  d4(ju>)4  +  d3(jco)3  +  d2(jw)2  +  dx(jw)  +  dQ 

The  I.  equation  was  programmed  and  solved  on  the  digital  computer. 
The  program  was  written  so  that  the  type  aircraft  could  be  changed,  i.e. 
all  stability  derivatives  were  retained  in  the  computation  of  xl  through 
x7  and  zl  through  z4.  The  gust  characteristic  length  (CL)  was  set  at 
the  beginning  of  the  program.  To  facilitate  handling  of  the  variables 
Kl,  K2,  and  K3  their  coefficients  were  stored  in  3  dimensional  arrays 
according  to  the  exponents  of  the  variables.  For  example,  the  coefficient 
of  K12K2K33  would  be  stored  in  an  array  at  location  (3,2,4) .  The  value 
stored  in  this  location  was  printed  out  as  PI (2,1,3)  =  (value  stored)  . 
The  arrays  were  manipulated  to  solve  for  the  performance  index  in  terms 
of  Kl,  K2,  and  K3. 

A  ndnimum  value  of  PI(K1,  K2,  K3)  was  determined  by  a  steepest 
descent  method.  A  starting  point  for  the  minimization  procedure  was 
obtained  by  computing  values  of  the  performance  index  as  Kl,  K2,  and  K3 
varied  over  a  large  range.  The  point  where  the  performance  index  was 
a  minimum  was  then  used  as  a  starting  point  for  the  steepest  descent 
method.  Several  other  starting  points  were  utilized  to  help  ensure  that 
a  proper  ininimum  was  reached. 
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II.  Flow  Chart 


VAK 


SET  CL,  Kl,  K2,  AND  K3 
ARG(l)  =  Kl?  ARG(2)  =  K2;  ARG(3) 


=3) 


/COMPUTE  AND  SET  VALUES  OF  d  AN&i 
V  c  TO  SOLVE  THE  INTEGRAL    J 


/SEE 


EE  PAGE  41  OF  THIS  APPENDIX 
MEANING  OF  d  AND  d 


D 


CALL  SUBROUTINES  TO  CALCULATE 
PERFORMANCE  INDEX,  SEE  SECTION  III 
OF  THIS  APPENDIX  FOR  SUBROUTINE 
EXPLANATION 


WRITE  OUT  PERFORMANCE  INDEX  NUMERATOR 
AND  DENOMINATOR  IN  THE  FORM: 
PI(I,J,K) 
WHERE  I, J,  AND  K  ARE  THE  EXPONENTS 
OF  Kl,  K2,  AND  K3  RESPECTIVALLY . 
1(2,  1,  3)  =  .06  WOULD  BE  .06Kl2K2K3l 


COMPUTE  THE  PARTIALS  OF  THE  PERFORMANCE 
INDEX  NUMERATOR  AND  DENOMINATOR 
WITH  RESPECT  TO  Kl,  K2,  AND  K3 


CALL  SUBROUTINE  FUNCT  WHICH  DETERMINES 
THE  VALUE  OF  THE  PERFORMANCE  INDEX 
AND  COMPUTES  THE  GRADIENTS  IN  THE 
Kl,  K2,  AND  K3  DIRECTIONS  FOR  A 
GIVEN  Kl,  K2,  AND  K3 .   THIS 
SUBROUTINE  USES  THE  PARTIALS 
COMPUTED  IMMEDIATELY  ABOVE. 
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PRINT  PERFORMANCE  INDEX, 
Kl,  K2,  K3,  AND  GRADIENTS 
IN  Kl,  K2,  AND  K3  DIRECTIONS. 


<5 


ALL  GRADIENTS  LESS  THAN  10~6 


< 


NO 


STEP  SIZE  LESS  THAN  10 


YES 


3uT\>- 


NO 


YSTEP  IN  DIRECTION  01 
V  STEEPEST  DESCENT 


<  ANYMORE  CC 
cL,  Klf  K2 


CONSTANTS 


AND  K3]^ 


NO 


YES 


YES 
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III.  Subroutines 

ZEPD3   -  Places  zeros  in  a  3x3x3  matrix. 

ZER04   -  Places  zeros  in  a  4x4x4  matrix. 

MULT1   -  Multiplies  a  first  order  polynomial  by  a  first  order 
polynomial  and  stores  answer  in  a  3x3x3  matrix. 

MJLT2   -  Determines  the  product  of  a  first  order  polynomial 
times  a  first  order  polynomial  times  2.0  and  stores 
answer  in  a  3x3x3  matrix. 

MULT3   -  Multiplies  a  first  order  polynomial  by  a  second  order 
polynomial  and  stores  answer  in  a  4x4x4  matrix. 

SUBT1   -  Subtracts  a  third  order  polynomial  from  a  third  order 
polynomial  and  stores  answer  in  a  4x4x4  matrix. 

SUBT2   -  Subtracts  a  second  order  polynomial  from  a  second 

order  polynomial  and  stores  answer  in  a  4x4x4  matrix. 

STUFEP  -  Determines  the  product  of  a  second  order  polynomial 
times  a  third  order  polynomial  times  a  constant  and 
stores  the  answer  in  a  6x6x6  matrix.  Separate  6x6x6 
matrices  hold  the  numerator  and  the  denominator  of 
the  performance  index. 

FUNCT   -  Takes  the  performance  index  numerator  and  denominator 
arrays,  the  specific  values  of  Kl,  K2,  and  K3,  and  the 
arrays  holding  the  partials  of  the  performance  index 
numerator  and  denominator  and  calculates  the  value  of 
the  performance  index  and  the  gradients  in  the  Kl,  K2 
and  K3  directions. 
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readily  measurable  normal  acceleration  factor  (n)  and  pitch  rate  (q)  can  be  used 
to  provide  a  and  a  information  indirectly. 
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